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Abstract. - Dielectric measurements on molecular liquids just above the glass transition 
indicate that alpha relaxation is characterized by a generic high-frequency loss varying as 
u -1 / 2 , whereas deviations from this come from one or more low-lying beta processes [Olsen et 
al, Phys. Rev. Lett. 86 (2001) 1271]. Assuming that long- wavelength fluctuations dominate 
the dynamics, a model for the dielectric alpha relaxation based on the simplest coupling between 
the density and dipole density fields is proposed here. The model, which is solved in second 
order perturbation theory in the Gaussian approximation, reproduces the generic features of 
alpha relaxation. 



Introduction. Alpha relaxation is the slowest and the dominant relaxation process 
in viscous liquids, closely linked to the glass transition because it sets the time scale for 
equilibration. Two of the most important unsolved problems regarding viscous liquids concern 
properties of the alpha process. One problem [1] is to explain the temperature dependence 
of the alpha relaxation time r, which for some liquids increases by more than a factor of 10 
when temperature is lowered by just 1%. Here we address the second problem: how can the 
observed shape of the alpha loss peak be explained? 

It was recognized already more than 40 years ago that the dielectric alpha process is always 
asymmetric, with almost Debye behavior at frequencies lower than the loss-peak frequency 
whereas there is a non-trivial power law loss, e"(w) cx tu~ n , above the loss peak [e(uj) = e'(uj) — 
ie"(ui) is the frequency-dependent dielectric constant, e" the dielectric loss]. In many cases n 
was found to be close to 1/2 and over time several theories were proposed predicting n = 1/2 
[2]. As measurements accumulated apparently contradicting any such universality, it became 
popular to represent data by the stretched exponential dipole time-autocorrelation function (oc 
exp[— (i/r)' 3 ] where t is time and (3 < 1 a parameter). Via the fluctuation-dissipation theorem 
this function reproduces the Debye behaviour at low frequencies, and it predicts an asymptotic 
power-law decay above the loss peak with n = (3. Often, however, the high-frequency power 
law applies only just above the loss peak frequency and a "wing" appears at higher frequencies. 
Nagel and co-workers proposed that, when scaled in an unusual way, the alpha dielectric 
loss follows a universal curve which includes the wing [3]. In 1997 at the 3rd International 
Discussion Meeting on Relaxations in Complex Systems it was suggested by Olsen, however, 

© EDP Sciences 



2 



EUROPHYSICS LETTERS 



that beta processes may play a role at much lower frequencies than previously expected and 
that this could explain the wing [4]. This was confirmed in experiments by Lunkenhcimer 
and co-workers and subsequently by other groups, showing that the wing indeed develops 
into a separate relaxation process after long-time annealing right below the glass transition 
[5]. A publication from 2001 giving data for ten molecular liquids [6] presented indications 
that n = 1/2 is a universal high-frequency exponent in the following sense: Equilibrium 
viscous liquids approach this behaviour as temperature is lowered and effects of beta processes 
gradually become negligible; in some cases, however, this may happen only when the alpha 
relaxation time is months or longer, making it impossible to verify or falsify the conjecture 
that n = 1/2 is universal. 

It is not generally accepted that n — 1/2 is a generic high-frequency exponent of alpha 
relaxation. Nevertheless, it makes sense to look for simple theories which can reproduce the 
proposed generic behaviour; the exponent n — 1/2 ought to be easier to explain than any other 
non-trivial exponent. For reasons not to be detailed here, the old models all have problems. In 
a recent paper it was suggested that a long-time-tail mechanism might be at the root of generic 
alpha relaxation [7]. Two immediate questions arise if this were correct: 1. The mechanism 
should apply a few decades above the alpha loss peak frequency, i.e., at times shorter than the 
alpha relaxation time; 2. the long-time tail of the relevant velocity autocorrelation function 
must be negative. The latter problem is solved by assuming stochastic dynamics, whereas 
a solution of the former is suggested by the solidity of viscous liquids [7]. Reference [7] did 
not discuss a specific model for the alpha relaxation. This is done below, where a model 
implementing the long-time-tail mechanism is proposed and solved for e(uj) in the simplest 
approximation. 

Model. - Consider N molecules in volume V. If n is the normalized molecular dipole 
vector, the basic degrees of freedom are taken to be the density field p(r) and the dipole- 
dcnsity field n(r). It is convenient to transform into k-space by defining the following sums 
over all molecules 



n k 



^5' 

3 



3 



N ■ 

3 



If [3 is the inverse temperature, the Hamiltonian (free energy) for a Gaussian model Hq is given 
by (3Hq = ^ k /9k/f-k/25(k) + X)k n k ' n -k/2e(k) where S(k) is the static structure factor 
and e(k) the dimensionless wavevector-dependent static dielectric constant. The simplest 
interaction term is H mt cx J n 2 (r)p(r)dr. We shall assume that the dynamics are dominated 
by the long-wavelength behaviour of the fields. Denoting the k — > (bare) limits of S(k) and 
e(k) by A and B respectively, in terms of a dimensionless coupling constant A the Hamiltonian 
becomes 

1311 = P*P- k + 5b S nk ' n k + 77^ E nk ' n k'P-k-k' • (2) 

k k v k,k' 

This defines an "ultra- local" field theory, i.e., one where equal-time fluctuations are uncor- 
rected in space. A more realistic Hamiltonian has gradient terms, but ignoring these is 
consistent with the assumption of long-wavelength dominance. Since p is always equal to 
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y/N, it is not a dynamic degree of freedom and all sums are to be understood excluding terms 
with po. 

Langevin equations define the dynamics: 



dp. 



n k = -r n (fc)^^ + r, k (i), (3) 

OTl_ k 

where £ k (/;) and ?7 k (t) are standard Gaussian white noise terms. Because density is conserved 
one expects r p (fc) cx k 2 . We define Y p (k)/A = Dk 2 . The dipole-density field, on the other 
hand, is not conserved. Consistent with long- wavelength dominance we identify r n (fc) with 
its k — * limit by writing T n (k)/B = V. Substituting eq. © into eq. J2J now leads to 



p k = -Dk 2 (p k + ^i^n k+k / n k , ] + £ k (t) 



n k 



N V 



r ( nk + W ? nk+k ''°- k ' ) + Vk (t) ■ (4) 



Equations (J2J and Q imply that to zero'th order in A the autocorrelation functions are given 

by 

(0k(O)p-k(<)) o = Ae- DkH 
(n k (0) • n_ k (t))„ = 35e- n . (5) 

We proceed to calculate the dipole autocorrelation function by applying the following theorem: 
If Qi obeys a Langevin equation of the form Qi = —fxdiH + one has (no sum over z) 

d 2 /cft 2 (Qi(0)Qi(t)) = p 2 (diH(0)diH(t)) [7,8]. To second order in A the dipole autocorrelation 
function in the Gaussian approximation thus obeys 

r- 2 ^(n k (0)-n_ k (t)) = <n k (0).n_ k (t)) 

+ ^E^'W^-wio^^^'io' (6) 

k' 

The sum is equal to 3ABVe~ rt (Dt)~ 3 / 2 /327T 5 / 2 . In terms of the dimensionless time t = Tt, if 
the constant 3A 2 AB 3 (T / D) 3 / 2 V/N8tt 5 / 2 is denoted by A, the dipole autocorrelation function 
$(t) = (n k (0) • n_k(f)) is the solution of 

$ = $ + Ae-*> 3/2 . (7) 
The general solution to this equation obeying <I>(i — * oo) = is 

$(t) = Ay sinh(? - t) e -*'F- 3/2 d? + C 2 e~*, (8) 
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which can be expressed in terms of the error function. In dimensionless units the dielectric 



This function has the required asymptotic behaviours: e"(u) cx ui for lot <C 1 and e"(oj) cx 
lj~ x / 2 for ujt 3> 1. There is one overall scaling parameter and one non-trivial parameter 
determining the loss-peak shape around the maximum. 

Discussion. - It has been shown that generic features of the alpha process are reproduced 
by a model based on a conserved scalar field and a non-conserved vector field with the simplest 
possible interaction term. The basic assumption of the model is that the dynamics are dom- 
inated by long-wavelength fluctuations. This justifies taking T p (k) oc k 2 and that S(k) and 
e(k) may be replaced by their k — > limits. Note that long-wavelength dominance is unrelated 
to possibly having a diverging length scale upon cooling. It simply means that the coherent 
diffusion constant D is much larger than estimated from r and the average intermolccular 
distance a (D » « 2 / T ) - whenever this applies, density fluctuations on length scales of order 
10a-100a take place on the alpha time scale and influence the dielectric loss, even above the 
loss peak frequency. The resulting dipolar dynamics have no length-scale dependence. 

An argument for long- wavelength dominance of the dynamics is given in a recent publica- 
tion [7] , where it is argued that the solidity of highly viscous liquids implies that the coherent 
diffusion constant D is much larger than the incoherent diffusion constant D s [defined from 
the long-time mean-square displacement of one molecule] . The reasoning, which assumes that 
the coherent diffusion constant is frequency independent, is briefly summarized as follows. At 
high viscosity molecular motion occurs via "inherent dynamics", i.e., jumps between potential 
energy minima [9]. Momentum conservation is irrelevant when relaxation times are in the 
second range or longer (the liquid container walls contribute any required momentum) . Thus 
if X is the sum of all molecular x-coordinates, it makes sense to define a "center of mass 
diffusion constant" by Dqm = lim^oo limy^oo (AA 2 (i)) /2Nt, where N/V — constant as 
V — > oo . For the high-frequency limit Z?cm(oo) one finds S(0) = Dcm{oo)/D which, if D is 
assumed to be frequency independent, implies £>cm(oo) <C D. Moreover, it is argued that 
D s (oo) ~ Dcm(°o) and D s (oo) 3> D s , where D s is the dc limit of D s (u). Naively one expects 
D s ~ a 2 /r, but the frequently observed breakdown of the Debye-Stokes-Einstein relation [10] 
implies that D s ^> a 2 /r. In either case, combining all these inequalities leads to O > a 2 /r. 
If one assumes that T p (k) cx k 2 applies when 1/T p is longer than the characteristic time 
t c = a 2 /D, we find that, since t c <C r, T p (k) cx k 2 applies for a range of rates larger than 
the loss peak frequency. - Note that, even if diffusion constant decoupling does not apply 
and D <~ D Sl long- wavelengths dominate the dynamics whenever the Debye-Stokes-Einstein 
relation is violated. 

The model only applies for extremely viscous liquids (with viscosity at least 10 10 times 
larger than that of ambient water), because only in this limit do liquids have the solidity 
property. Experimentally, this is also where the w -1 / 2 high-frequency behaviour is asymptot- 
ically approached. We conjecture that in the extremely high viscosity limit (but still in the 
equilibrium liquid phase) an uj~ x I 2 high-frequency universality always appears. This, in fact, 
should apply to any linear response function. 

A few final comments: 

1. The assumption that T p (k) cx k 2 cannot - despite its hydrodynamic flavour - apply all 
the way down to k = 0. At very large wavelengths a k- independent density "decay channel" 




Ay/ir and 
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must appear, because otherwise there is no bulk volume relaxation on the alpha time scale 
and the glass would have the same compressibility as the equilibrium liquid. Thus a more 
accurate model would assume r p (/c) = To+Dk 2 (which, however, only slightly perturbs model 
predictions if Tq ^ 1/r as required by experiment). 

2. The Hamiltonian of eq. J5J is "bottomless" in the sense that it allows arbitrarily low 
energies [11], implying that the canonical distribution is not normalizable. This inconsistency 
is readily resolved by assuming that there are higher even order terms in n(r). If these terms 
are small, they barely influence the equilibrium autocorrelation functions, but still ensure that 
the fields do not eventually run off to infinity. 

3. According to the model the generic uj^ 1 / 2 high-frequency decay of the alpha process is 
specific to three dimensions. In four dimensions, for instance, the sum in eq. (0 varies with 
time as t~ 2 , leading to an autocorrelation function which varies logarithmically with time as 
t -> 0. 
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APPENDIX 

In order to calculate the Laplace transform of — $, we first note that the second term of eq. 
JHJ leads to the C2-term of eq. The first term requires calculating the following integral 
(where s = icu and tildes are left out for simplicity of notation): 



o 



cosh(i' - t)e~ 1 ' tl^^dt' = I dt 
We split I into three terms: 



°° e t'-t + e *-f 



i pOO pOO i p OO pOO 

- / dt e-< s+1 >* / t'- 3 / 2 dt' + - dt e- {s -^ / t'- 3 / 2 dt' + (*), 



2 7o 



where 



Thus 



-| poo poo 

(*) = ~ / dt e^- 1 ^ / (e" 24 ' - l)i'~ 3 / 2 ^'. 

^ .In .1 + 



1= / e- {s+1 H- 1 / 2 dt+ / e-^H-^dt+i*). 
Jo Jo 

Since e^ t t^ 1 l 2 dt = y/ir we get 

/ = v^(i + s)- l/2 + v^(s - iy 1/2 + (*). 



The last term is calculated as follows: 



(*) 



1 - s 



( S -l)t / ( e _2f _ 1 ) t /-3/2 dt / 



t=0 



1 1 

2 1 - s 



-(-D*(X_ e'^r^dt. 



Thus 
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or 



where 



(*) = \^— g QT(l - e- 2t )r 3 / 2 dt - £ (e-^- 1 )* - e -< s+1 )*) t-*l 2 dt 
(*) = T l_(F(0,2)-F(*-l,a + l)) 



1 f 00 

F(o,6) = i / (e- at - e" 6 ') r^ 2 dt. 

2 Jo 

Calculating this function is straightforward: 

F{a,b) = (e- at - e- bt )(~t-^ 2 ) - (ae- at -be- bt )t-^ 2 dt = (P{b)-^a), 
L J t=o J a 



where 



4>{x) = / xe- xt r 1/2 dt = vW^- 
Jo 

Summarizing, 

-4 = (i + s)- 1 ' 2 + ( fl - + _L_ (V2 - Vo - + V^T 

y7T 1 - s V 



or 



s 1 



i-s vm 1 v2 + ^/TTs 
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